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The main goal of ultrarelativistic heavy ion collisions (URHIC) is to create strongly 
interacting condensed elementary particle matter, quark-gluon plasma (QGP), and to study 
thermodynamics of this new phase of matter. The existence of the QGP phase has been 
proved by ab initio calculations in lattice QCD |TJ. 

High energy particle physics usually aims at producing as few particles as possible in 
order to create a simple system for the measurements. The URHIC also aim at a simple 
system but by producing as many QCD quanta as possible. In hadronic collisions the 
produced final state is dilute, further interactions negligible, and the system essentially 
freely streaming. In URHIC, one would ideally produce a system with a maximal amount of 
final state interactions, i.e. a large enough thermal parton system with negligible mean free 
paths. 

From the point of view of QGP formation, as high cms-energies yfs and as large nuclei 
A as possible are preferable. So far, the most energetic nucleus-nucleus collisions have been 
the Au-Au collisions at y/s = 5 AGeV at the AGS0 in the Brookhaven National Laboratory 
(BNL), and Pb-Pb collisions at yfs = 17.6 AGeV at the SPSQin CERN. A very strong J/tf- 
suppression observed in the central Pb-Pb collisions at the SPS by the NA50 collaboration 
suggests that the QGP phase may have already been reached at the SPS, and, more 
importantly, observed. 

In future heavy ion collisions, Au-Au at y/s = 200 AGeV at the BNL RHICQ, and Pb-Pb 
at y/s = 5500 AGeV at the CERN LHC/ALICEfl, initial energy densities far beyond the 
critical densities of QGP formation will be produced. It is also expected that the system 
in the central rapidity region of the RHIC and LHC nuclear collisions will stay in the QGP 
phase for several fm/c, thus giving better possibilities for observing characteristic signals of 
the QGP like J/\l/-suppression, thermal production of dileptons and photons, strangeness 
enhancement, collective flow, disoriented chiral condensates, etc. For a review, see Ref. @. 

Particle and transverse energy production in the central rapidity region of URHIC can 
be treated as a combination of perturbative (hard and semihard) parton production and 
non-perturbative (soft) particle production. By "hard processes" one usually means clearly 
perturbative processes with momentum or mass scales of the order of several or tens of GeV, 
while "semihard" here refers to QCD-processes where partons with transverse momenta of 
a few GeV, "minijets," are produced. In this article I will discuss the role of perturbative 
parton production as the early initial conditions for the QGP j|, [5], @] i n URHIC at LHC 
and RHIC. In particular, some recent results of initial parton production [^, 0] are reviewed 
in Sees. 2 and 3. I will also briefly discuss other approaches [§, [1(J in Sec. 4, and further 
evolution of the QGP in Sec. 5. 

By definition, the semihard processes lie at the border-line of hard and soft physics. 
Real hadronic jets have been observed in pp collisions from p^ > 5 GeV [JXTj] up to p^ ~ 



440 GeV [12]. The minijets with px ~ 1-..2 GeV are part of the underlying event rather 
than distinguishable jetsQ. In nuclear collisions, where thousands (hundreds) of minijets are 
expected to be produced in the central rapidity unit at the LHC (RHIC) j|, ||, detection 

2 AGS = Alternating Gradient Synchrotron 
3 SPS = Super Proton Synchrotron 
4 RHIC = Rclativistic Heavy Ion Collider 

5 LHC = Large Hadron Collider; ALICE = A Large Ion Collider Experiment 

6 In this sense "minijet" is not a very good name for a semihard parton; it is not a "jet" at all. 
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of individual semihard partons becomes impossible. 

Even though the semihard partons are not observed as jets, their production should 
fall well within the scope of perturbative QCD (pQCD) since p T 3> A QCD ~ 200 MeV. 
Ultimately, to how low values of p^ one can push the validity of pQCD, is a question of 
convergence of the QCD perturbation series. The recent results from HERAfj at DESY 
TJ, |nj] indicate that the behaviour of parton densities can be explained within pQCD down 
to the low scales Q 2 ~ 1 GeV 2 . This supports applicability of pQCD to semihard parton 
production as well. 

In the division into semihard and soft particle production, the key feature is that one 
is able to compute the semihard parton production from pQCD, given that the parton dis- 
tributions of colliding hadrons and nuclei are known. On top of this, the non-perturbative 
particle production can be modelled e.g. through strings [|l|, [Tj| [17|, |18| , or in URHIC 



perhaps also through a decaying strong background colour field 1T9L |20[ . With increasing en- 



ergies however, the semihard QCD-processes are expected to become increasingly important, 
particularly in URHIC. This is due to the following reasons: 



With increasing cms-energies, the events tend to become more "jetty" or, rather, more 
"minijetty" , and in pp and pp collisions the rapid rise of the total and inelastic cross 
sections with the \fs can be explained by the copious production of semihard partons 



with transverse momenta p^ > po ~ 1...2 GeV [£1], Perturbative parton production 
in A A collisions scales as ~ A 4 / 3 j|, |6j, while the soft component scales more like ~ A, 
so for large nuclei the importance of semihard partons should be further increased. 

In the deep inelastic ep scatterings (DIS) at HERA it is observed that the structure 
function F%(x, Q 2 ) has a steep rise at small values of Bjorken x, at x ^ 0.01, persisting 
down to scales Q ~ 1 GeV [^3[ [14], [13| . The quark-antiquark sea is generated by 



emission from the gluons, so the gluon distributions have this rapid rise as well. In 
minijet production in the central rapidities, the dominant gluonic processes at ~ 
2 GeV probe the parton distributions typically at fractional momenta x ~ 2p^ / yfs ~ 
7 x 10~ 4 in the LHC nuclear collisions. For the RHIC, where typically x ~ 0.01. ..0.02, 
the rise of the parton distributions will not cause such a big effect at the central 
rapidities. 

Time scale for producing partons and transverse energy into the central rapidity region 
by semihard collisions is short, typically 7"h ~ l/po ~ 0.1 fm/c, where po ~ 2 GeV is 
the smallest transverse momentum included in the computation. The soft processes 
are completed at later stages of the collision, at r s ~ 1/Aqcd ~ 1 fm/c. If the density 
of partons produced during the hard and semihard stages of the heavy ion collision 
becomes high enough - there are indications that it will - fusions start to occur, and 
a saturation in the initial parton production can take place in the perturbative region 
[0, H. As a result, softer particle production will be screened. The fortunate 



consequence of this is that a larger part of transverse energy production in the central 
rapidities is computable from pQCD at higher energies and the relative contribution 



7 HERA = Hadron Electron Ring Anlage 
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from soft collisions with pt ~ 2 GeV becomes smaller. Typically, the expectation is 
that at the SPS the soft component clearly dominates, and at the LHC the semihard 



component is the dominant one || |18[ . At the RHIC both components should be 
taken into account. 

The importance of copious semihard parton production in an early formation of the 
QGP in URHIC was first addressed almost 10 years ago, by Blaizot and Mueller [H, and by 
Kajantie, Lindfors and Landshoff [[|. In terms of eikonal models for hadronic interactions, 
minijet production had already been emphasized some years before, in Refs. j2l| , to explain 



the rise in the total and inelastic cross sections of pp cross sections beyond the CERN 
ISR energy range. Models including both strings and semihard parton production were 
introduced around the same time . The work || was later improved by including partonic 
cross sections in lowest order (LO) pQCD |J and with a better treatment of the rapidity 
acceptance. Consequences of the HERA results for minijet production in URHIC were 
studied in 0. Minijet production in AA collisions has also been studied by Calucci and 



Treleani from the beginning of the 90's |p7[| , and later by Xiong and Shuryak |Pq| . In the 
end of the 80's, systematic efforts began to construct event-generators for URHIC, which 
would be based on the increasingly important pQCD component. As a result, HIJING by 



X.-N. Wang and Gyulassy and, Parton Cascade Model by Geiger and Miiller [^9] were 
constructed. Also other event generators have options for simulating the URHIC at the 
RHIC and LHC energies §§,§. 

Treatment of semihard Q CD-scatterings is usually based on collinear factorization, where 
the hard parton-parton scatterings are factorized from the universal parton distributions 
at a perturbative scale ~ > Aqcd- The parton distributions then contain the non- 
perturbative experimental input. A novel approach to semihard parton production, not 
based on the collinear factorization nor independent scatterings, but on a consideration 



of a classical gluon field, has been suggested by McLerran and Venugopalan in 1994 [|0 
Not based on collinear factorization either, and perhaps directly related to the gluon field 
approach, minijet production in a BFKL approach was studied recently ||. These topics 
will be discussed in Sec. 4. 



1 Minijet production in AA collisions 

The idea of multiple production of semihard gluons and quarks in pp and especially in AA 
collisions is based on a picture of independent binary parton-parton collisions. The key 
quantity is the integrated jet cross section: 

<7jet(Vs,Po,Ay) = - dpldy 1 dy 2 Yl dy 2 x 1 f i/N (x 1 ,Q) x 2 f j/N (x 2 ,Q)- I (s,t,u), (1) 
l Jp 2 Ay ij k i = J at 

q,q,g 

where x\ t2 are fractional momenta of the incoming partons i and j. Parton- level quantities are 
indicated by hats, and fi/N^x, Q) are the parton distributions in N (= p, A) at a factorization 
scale, chosen as Q = p^. The normalization factor 2 comes from the fact that, in the lowest 
order (LO) pQCD, there are two partons produced in each semihard subcollision. In the 
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eikonal models for pp collisions |22| the ratio <Tj Ct / 0"mciastic can be interpreted as the average 



number of semihard collisions in one inelastic event. The explicit numbers I will be quoting 



in the following || are obtained with the GRV-LO parton distributions [[3T]] exhibiting a 
small-x rise similar to that in the HERA data. More detailed formulation can be found in 
Refs. 0, ||, and a numerical evaluation of Eq. (p]) with other parton distributions in Ref. 

ED- 

The above formula is defined in the lowest order, da/dt ~ a 2 . Often a constant factor 
K ~ 2 is used to simulate the effects of next-to-leading-order (NLO) terms. Studies of the 
NLO jet cross section da / '(dp^dy) |32| show that (with a scale choice Q = Pt and with a 



jet size R ~ 1) this may be a reasonable approximation |p3 |. Strictly speaking, however, 
a theoretical .fT-factor can only be defined for quantities where a well-defined, infrared- 
safe measurement function can be applied |32]. For ^-production in nuclear collisions, 



an acceptance window in the whole central rapidity unit defines such a function, but for 
this acceptance criterion, and for p^ ~ 2 GeV, the exact NLO contribution has not been 
considered yet. For consistency reasons, however, there is no K-factor included in the explicit 
results I will discuss here. 

A first estimate of the average number of produced semihard partons in a rapidity window 
Ay with pt > po in an AA collision at a fixed impact parameter b can be obtained as [|j] 

NAA(b,y/s,p Q ,Ay) = 2T AA (h)a }et (y/s,p ,Ay), (2) 

and the average transverse energy carried by these partons as [|J 

E^ A (b,^,Po,Ay) = T AA {b)a iet {^,p )(E T ) Ay , (3) 

where T AA (b) is the nuclear overlap function. The normalization is / d 2 hT AA (h) = A 2 , and 
since T AA ~ A 4 / 3 , it describes the typical scaling of hard processes in nuclear collisions. 
For large nuclei with Woods-Saxon nuclear densities, T AA (0) « A 2 /(nR A ). The acceptance 
criterion Ay will be \y\ < 0.5, and corresponding cuts will be made for yi and ?/2 @- In 
Eqs. (H) and (||) above, T AA {h)a- iet is the average number of semihard collisions and (Et)a v 
is the average transverse energy carried by the partons produced in each of these collisions 
into Ay. Parton-flavour decomposition and rapidity distributions can be found in ||. 

In Figs. [1], the integrated jet cross sections and the first £/r-nioments for y/s = 5500 
and 200 GeV are shown as functions of the smallest transverse momentum po included in 
the computation. These quantities naturally depend strongly on the choice of po since the 
subcross sections diverge as da/dt ~ p^ 4 when p T — > 0. Since p is a parameter that decides 
the division between soft and hard physics, by definition some phenomenology is needed to 
fix its value. One way to have control over po is to study an eikonal model [p^j , where the 



lower limit of po (and simultaneously a possible .fT-factor) is controlled by the rise of the cross 
sections. In addition, by fitting measured charged particle spectra in hadronic collisions at 
high energies, it is possible to extract a value for p in connection with a string model [^6| . 
Also, convoluting the partonic cross sections with fragmentation functions of each parton 



flavour into charged pions and kaons [35[] gives a simultaneous handle (in the LO) on the 
.fT-factor at large pt and on p at the small pt part of the charged particle |>r-spectra. These 
procedures result in values p ~ 2 GeV. 
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Figure 1: The integrated minijet cross section aj e t(po, y/s, \y\ < 0.5) and the first .Ex-moment 
c r ict(Vs,po){E T ) A y as functions of p at \y\ < 0.5 at LHC (y/s = 5500 GeV) and RHIC (y/s = 200 
GeV) from Eqs. @ and @. The dashed curves are the gluon contributions, the dot-dashed curve 
shows the transverse saturation limit for A = 208. Shadowing is not included and K = 1. 

In URHIC then, if sufficiently many partons (mainly gluons) are produced in the very 
beginning of the collision, the partons (within Ay) start to overpopulate the available nu- 
clear transverse area ~ kR 2 a , and final-state fusions become important || [24j| . In this 
case, further production of partons, especially the transverse energy production, becomes 
screened. Let us assign an effective transverse area (uncertainty) 7r/p T for each parton 
produced within \y\ < 0.5. Since the partons dominantly have ~ p , we can estimate 
0] that a saturation in the semihard parton production in URHIC should happen when 
N AA (p ,y/s, \y\ < 0.5)-k/pI>-kR 2 a , i.e. when 

^(V~s,Po, \y\ < 0.5) > ~ A~ 2/3 pl, (4) 

so the larger the nucleus is, the earlier the saturation will occur. The saturation value is 
plotted in Fig.|l| for A = 208. The saturation of the minijet production cross section can be 
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expected near p^ ~ 2 GeV for the LHC and near ~ 1 GeV for the RHIC. Therefore, for 
the LHC, the bulk of transverse energy in the central rapidity unit is expected to come from 
the semihard processes alone, so po = 2 GeV is a reasonable choice for the LHC. The eikonal 
model and the fragmentation function analysis indicate that a choice po = 1 GeV would 
result in an overestimate of the (inelastic) cross sections and the charged particle pp-spectra. 
Instead of trying to fit this value any better for the RHIC, we will consider here the initial 
conditions for QGP production with the same value po = 2 GeV as for the LHC. 

Self-consistent screening of initial parton production can also be modelled in terms of 
a dynamical, medium-induced screening mass p5|| . The midrapidity partons are produced 
at times r ~ 1/pt, so the large-pr partons are produced first, and partons with smaller p-r 
are produced later. In a medium of high-pr partons, a dynamical screening mass (electric, 
static) m g []34]| is generated, and this mass then screens the production of partons with 
smaller p^. If the density of produced partons is high enough, i.e. if the -y/s and A are large 
enough, the screening mass grows fast enough, causing a saturation of the parton cross section 
already in the perturbative regime pt 3> A QCD . In practice, since the transverse part of the 
screening is not known, we have modelled the screening by simply making the replacement 
t{u) — > t(u) — m 2 in the partonic cross sections. Saturation in this approach coincides with 
Eq. (H), and, by computing the first moment of the px-distributions, one also concretely 
observes how the bulk of transverse-energy production is obtained from p^ >2 GeV at the 
LHC (see J3§). 

To conclude this section, I will discuss nuclear parton distributions. In the compu- 
tation presented above, these are approximated as fi/A(x,Q 2 ) = Afi/ p (x,Q 2 ). Clearly, 
this gives a fair first estimate. It is, however, experimentally known that there are nu- 
clear effects to the parton distributions. In the DIS measurements |3B], [J7|, the ra- 
tio F^(x, Q 2 )/F 2 D (x, Q 2 ) has been measured, and four regions in x can be distinguished: 
depletion at x<0.1 called "nuclear shadowing", enhancement at 0.1<x<0.3 called "anti- 
shadowing", depletion at 0.3^ x ^0.7 called "emc-effect" , and cumulative enhancement at 
x — > 1 (and beyond) called "Fermi motion" . See Ref. |$9| for an extensive review of the 
experimental data and the various theoretical models. 

For the mid-rapidity minijet production, practically only the shadowing region is relevant. 
It is experimentally observed that the F 2 -ratio does not strongly depend on the virtuality 
Q 2 , so a scale- independent ratio may give a first estimate of the nuclear effects to the nuclear 
quark and antiquark distributions. However, as shown in [^T[] , it is not necessarily so for the 
gluons, but the QCD scale evolution |44] should be taken into account in more detail. In 
Fig. |2|, I show the evolution of ratios Rf(x,Q 2 ) = fi/A/Afi/ p separately for gluon, valence- 
and sea-quark distributions, and for F 2 , from |4T|. In this analysis, charge and momentum 
sum rules are incorporated with the DIS data p§| . 



and the dilepton data [40|. A reanalysis 
It is also becoming possible to 



with more modern parton distributions is being prepared [52 



extract the gluon distributions from the logarithmic derivatives of F 2 1 j 'F 2 



so further 



constraints for the badly known nuclear gluon distributions are available. The scale evolution 
of nuclear parton densities has also been studied by other people, see Refs. ||47|| . Finally, the 
role of the GLRMQ correction terms [24~, (which were included in [[H|]) to the DGLAP 
equations [44] should also be considered in more detail, perhaps along the lines in Ref. 



46 



but in the light of the most recent HERA-data 



To get an idea of the magnitude of 
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DGLAP+GLRMQ-evolution 




0.5 
0.001 



Figure 2: The ratios Rf(x,Q 2 ) = h, A (x,Q 2 )/Af l/p (x,Q 2 ) of i = g, V(= u v + d v ), S(= EIqs + 
Qs]), F% for A = 196, as functions of x at scales Q = 2,3, 5, 7, 10 GeV f4l| . The evolution is lowest 
order DGLAP, modified with the GLRMQ terms, and the data shown are from [40, 37] 
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\y\ < 0.5 


total 


9 


Q 


q 


LHC: 


ApbPb 


4740 


4349 


199.9 


191.2 




pPbPb 

_G/rp 


14170 


12960 


620.8 


588.8 


RHIC: 


^PbPb 


120.6 


99.58 


12.96 


8.102 




pPbPb 

£L/r£ 


320.7 


262.8 


36.06 


21.87 



Table 1: The average numbers and transverse energies of semihard partons at r = 0.1 fm/c with 
\y\ < 0.5 and > 2 GeV in central Pb-Pb collisions, as given by Eqs. (0) and (S). No shadowing 



nor K-factor is included and the GRV-LO |31j parton distributions are used. 



the shadowing effects in minijet production at the LHC and RHIC, I refer the reader to the 
computation in ET . 



2 Initial conditions for QGP at r = 0.1 fm/c from pQCD 

As described in the previous section, we obtain a first estimate of the initial conditions of 
QGP production by fixing the minimum as po = 2 GeV. In this way, we describe the initial 
conditions at r ~ l/p = 0.1 fm/c. The predictions for the average numbers and transverse 
energies of partons produced into the central rapidity unit in central Pb-Pb collisions at 
RHIC and LHC energies are summarized in Table 1 (with even too many decimals). Also 
individual contributions from gluons, quarks and antiquarks are indicated [^]. There are five 
quite straightforward but important observations: 



• Gluons strongly dominate the perturbative parton and transverse energy production: 
the initial parton system at r = 0.1 fm/c is about 90% glue in the LHC and 80% 
in the RHIC, so that at early times the QGP is actually gluon plasma to the first 
approximation. 

• As discussed in Sec. 2, the parton system at the LHC is transversally saturated at 
r = 0.1 fm/c, within the central rapidity unit. This means that production of more 
(softer) partons increases mainly the rate of final state fusions, not so much the trans- 
verse energy production. At the RHIC the same will happen but at a somewhat later 
stage which may not be entirely controllable by means of pQCD. In any case, it is 
demonstrated that the parameter po acquires a dynamical significance in the URHIC. 

• By assigning baryon number 1/3 (—1/3) for each quark (antiquark) produced, and by 
approximating the volume of the parton system by V = irRp h Ay/pQ = 13.4 fm 3 for Pb- 
Pb, the initial net baryon number density at r = 0.1 fm/c becomes n B _g = 0.21 fm -3 
for the LHC, and, 0.12 fm" 3 for the RHIC. Thus, even at these ultrarelativistic energies 
the initial net baryon number density is comparable to the nuclear matter density, 
0.17fm -3 , even beyond it at the LHC.0 Of course, the high initial net baryon density 

8 Note that by the time r = 0.1 fm/c, the Lorentz-contracted nuclear disks (the hard parts) are already 
receding; the transit time of the nuclei is tt ~ 2Ra/^ ~ 5 x 10 -3 fm/c for the LHC and 0.1 fm/c for the 
RHIC. 
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will dilute quite fast, but it is interesting that the pQCD computation gives such a 
high initial density. One would perhaps expect more baryon stopping at RHIC than at 
LHC. I should emphasize that the number presented here is only the perturbative part 
of net baryon number and that the increase for the LHC is entirely a gluon distribution 
effect because the perturbatively produced net quark number in the central rapidity 
region is mainly due to valence quark-gluon scattering. The typical values of x probed 
at the LHC are smaller than at the RHIC, so, because of the rapid rise of the gluon 
densities, there are many more gluons for the valence quarks to scatter with at the 
LHC than at the RHIC. 

To what extent can the initial system be considered thermalized? One immediately 
observes that there are far too few quarks and antiquarks as compared to gluons, 
so that the system cannot be in chemical equilibrium. How about the gluons alone 
then? In an ideal thermal gas of massless bosons the energy/particle is determined by 
€ g /n g = 2.7T eq . From the numbers for the LHC in Table |l]we can determine that 

j?g e pQCD 

^ = ^ocd =2.98GeV. (5) 



N 9 ^PQCD 
iv PbPb n 9 

In an ideal gas of massless gluons in complete thermal equilibrium, the temperature 
can be computed from e^ deal = 37r 2 /90 • 16T e 4 q , and for an ideal gas with an energy 
density ejf eal = eP QCD we find T eq = 1.10 GeV. We see that 

gPQCD gideal 

-I ~ -I — 9 7T (6) 

n g n g 

So, as far as the energy per particle is concerned, the gluon system is "thermalized" 
from the beginning at the LHC. At the RHIC, this is likely to happen somewhat later. 
One should, however, keep in mind that we did not consider isotropization at all here, 
and, that some uncertainty is connected with the assumption on an isotropization 
time. In order to estimate this time, a more detailed space-time picture of initial 
parton production is needed. Such a modelling is presented e.g. in Refs. |48], p9| . 

Finally, let us consider the initial net baryon-to-entropy ratio of the early QGP. For 
a thermal boson gas s = 3.60n, where n is the number density of thermal gluons, the 
total initial entropy (glue only) is S = 15900 for the LHC; so initially, at r ~ 0.1 fm/c, 
the net baryon-to-entropy ratio is (B — B)/S ~ 1/5000. For the RHIC it will be 
larger, about 2/1000. In the further evolution of the QGP, this number will increase 
somewhat due to the non-perturbative net baryon production; we estimated the final 
(B-B)/S ~ 8 x 10- 4 (9 x 10~ 3 ) for the LHC (RHIC) §. Even though we are studying 
the same phase transition as took place in the early Universe, we are still relatively 
far away from those extreme conditions regarding (B — B)/ S; there, the inverse of the 
specific entropy is ~ 10 -9 . 
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3 Beyond the factorized minijet approach 



The simple approach of Sec. 2 without shadowing effects in the parton distributions will 
fail at large rapidities when sufficiently high cms-energies and large nuclei are involved (see 
||). Minijet production at large y is a consequence of a few large-x partons scattering 
against a large number of small-x partons. Due to the rapid rise of the gluon distributions 
at small values of x, and without shadowing, the parton luminosity is high and the cross 
sections at large y are not suppressed enough for the assumption of independent scatterings 
to hold. Eventually, since Taa ~ A 4 / 3 , one will run into trouble in conserving energy and 
baryon number (which scale as ~ A) for a large enough A and y/s. In this case, coherence 
effects should be considered. However, within the central rapidity unit, where our focus 
is, the assumption of independent parton scatterings still works quite well: only energy 
~ 0(Naa x 2po) "C Ay/s is consumed, and only a fraction of the available number of 
partons will scatter. Naturally, the latter criterion also serves as a further constraint for the 
value of po in URHIC. I should also emphasize that the approach will work even better once 
nuclear (gluon) shadowing is included. 

The question of the validity of factorization in minijet production in URHIC is, however, 
an important one, and additional production mechanisms should be studied. During the 
recent years, quite a different approach to minijet production has been developed, not based 



on factorization. McLerran and Venugopalan |10[ have suggested a model where a large 
colour charge density of a large nucleus, travelling along the light-cone, generates a gluon 
field that is effectively classical. The idea then is that the gluon distribution function of the 
nucleus (perhaps even hadrons) could be computed in the region Aq CD <C k T <C fi, where 
is the transverse momentum of gluons, and the scale fi 2 is the area-density of gluons per unit 
rapidity. Connection of this approach to evolution equations, especially to the BFKL [ ^0|| , 
has been considered in ||54j| . The actual collision of two nuclei has been formulated in [ pSfl , 



and connection of the classical approach to Feynman diagrams has also been recently studied 



1 55], Also the connection to the BFKL- type minijet formula |52|, |9] is under investigation 
57j . It remains to be seen what the predictions of the gluon field approach eventually are in 
actual numbers for URHIC, and how well the gluon distributions and their QCD evolution 
in nuclei and hadrons are accounted for. In my opinion, it is very important to pursue this 
work into the direction outlined in [57|, where the applicability of the model to URHIC 
at the LHC and RHIC was studied in more detail. Ultimately, one could hope that while 
collinearly factorized minijets dominate parton production at pr ~ Po — 2 GeV, the region of 
Aqcd ~ Pt ~ Po (where also non-linear effects will eventually become important) could become 
better under control in terms of the BFKL and/or gluon field approach. 

Minijets can also be emitted from a BFKL ladder [|52][]. By assuming that the small- 
er increase in F% is entirely due to the BFKL physics ||50|| , the maximum transverse energy 
deposit in the central rapidity region due to the minijets from a BFKL ladder can be studied. 
This was done in ||. Since the increase in F 2 takes place only at x<0.01, the BFKL 
mechanism is not expected to be important in the RHIC nuclear collisions. 

Minijet production without high-pp tagging jets requires an introduction of unintegrated 

9 See the useful lecture notes by Del Duca for a derivation of the basic concepts and for original 
references. 
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parton densities, which evolve according to the BFKL equation |50[ (assumed to be homo- 
geneous |53|), and which can be normalized to the "known" gluon distribution xg(x,Q). 
Effectively, minijet production from a BFKL ladder can be considered as an 0(a s ), 2 — ► 1 
process where two virtual gluons fuse and form the minijet in the mid-rapidity. Techni- 
cally, this minijet is just one fixed rung of the BFKL ladder (the emitted gluons of which 
are strongly ordered in rapidity). In |J, we were unable to fix the overall normalization of 
the process, in the absence of any perturbative Born-level process to compare directly with. 
By comparing the BFKL computation with the factorized jet-production cross section (in 
the lowest order) at high pr, we were, however, able to argue that the transverse energy 
production from the BFKL minijets is still subleading at the LHC energies. 

Since the HERA results for the increase of F 2 at small x can be explained by the leading 
log(Q 2 ) and/or the leading log(Q 2 ) log(l/x) |49| approximations, the leading log(l/x) 
BFKL contribution is obviously not the dominant mechanism at the present values of x. 
Thus, my conclusion at this point is that the BFKL minijets certainly bridge the way to- 
wards softer physics at pr < Po ~ 2 GeV, but the initial conditions relevant for the early 
QGP formation in the LHC nuclear collisions seem to be dominantly given by the minijets 
computed in collinear factorization. 



4 Further evolution of the minijet plasma 

It is very important to understand when and how well the early parton system can be 
described in terms of hydrodynamics |59|]. Not only the thermal signals, but also global 
variables such as final transverse energy and total multiplicity will strongly depend on the 
onset of pressure and flow effects. Normally, in hydrodynamical calculations the measured 
final-state hadron spectra are used for getting an estimate of the initial conditions (see e.g. 
|58|1). These then, by definition, depend on the equation of state used. Semihard parton 



production could now provide the hydro codes with additional and independent information 
on the early initial conditions. Of course, validity of applying hydrodynamics to the early 
QGP ultimately depends on how completely and how fast the initial parton system reaches 
a state of local thermodynamic equilibrium. 

I have explicitly shown by using pQCD, how thousands (hundreds) of partons, mainly 
gluons, will be produced within r ~ 0.1 fm/c into the central rapidity unit of URHIC at 
the LHC (RHIC). For the LHC, the perturbative gluon system is shown to be transversally 
saturated and thermalized in the energy/particle sense. This indicates that already the early 
evolution of the QGP in the LHC could perhaps be described in a simple hydrodynamical 
approach ||59|| . Initially produced perturbative partons do not necessarily, however, have 
boost-invariant rapidity distributions ||. By making an extreme assumption of having a 
locally thermal (gluon) system at r ~ 0.1 fm/c, and by using the rapidity distributions as 
computed in || for the initial conditions, it is possible to study the deviations from a boost- 
invariant Bjorken picture. In |K| it is shown that the arising pressure gradients generate 
somewhat faster cooling of the (Q)GP as in the Bjorken picture. 

At the RHIC also the soft component in particle and transverse energy production should 
be taken into account. A possible modelling for this as a source term in Bjorken hydrody- 



namics, on top of the minijet initial conditions, was suggested in [20], and also effects of 
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dissipation and colour conductivity were studied. More complicated 3- dimensional systems 
with minijet initial conditions and density fluctuations have been considered in pll . 

Related both to the initial conditions and to the further evolution of the QGP, a more 
difficult and still open question is thermalization of quarks, i.e. how fast, if at all, quarks and 
antiquarks come into chemical equilibrium with gluons Phenomenologically, this subject 
is vital for the thermal electromagnetic signals |k|. Theoretically, the issue is space-time 
dependent phenomena in gauge theories, a virtually uncharted territory. 



5 Conclusions 

So far, the diverging minijet cross sections of Sec. 2 can be regulated in URHIC either 
by a dynamically determined cut-off parameter p , or by modelling in a screening mass. 
These are plausible but phenomenological considerations. Unless more theoretical progress 
in understanding pQCD parton scattering in the few-GeV range is made, these will remain so. 
A quantitative study of the NLO terms in the minijet cross sections, especially in transverse 
energy production, will certainly give us a better handle on the validity of pQCD at these 
scales. A resummation, as in the Drell-Yan dilepton case [|64 , will perhaps become possible 



eventually. As discussed in Sec. 2, the nuclear gluon distributions are also poorly known; 
ultimately they should be measured better. Coherent scattering and higher-twist effects 
should be studied in more detail, at large rapidities especially. 

The main results given in Sees. 2 and 3 are obtained by using collinear factorization and 
independent parton-parton scatterings. Minijets from the classical gluon field approach and 
from the BFKL ladder discussed in Sec. 4 are not based on factorization. By studying the 
upper limit of minijet production from the BFKL ladder, we were able to argue |J that the 
BFKL-mechanism should still be subleading in E^-production at the central rapidities, at 
least up to the LHC energies. Connection of the gluon field and BFKL approaches with the 
collinearly factorized minijet production should still be understood better [|57} . 

To conclude, it is clear that the analyses of semihard parton production discussed in 
this article can, and will be, sharpened in various ways. I have, however, no doubt that 
the semihard partons will play a major role in the formation and evolution of the QGP in 
central rapidities during the first fractions of fm/c of the LHC and RHIC nuclear collisions. 
Understanding the primary production mechanisms of partons will be the key to predicting 
and explaining the signals of the QGP and more global variables measured in the future 
heavy ion collisions at LHC and RHIC. 
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